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Abstract 

We present a simple, but efficient, way to calculate connection matrices between sets of indepen- 
dent local solutions, defined at two neighboring singular points, of Fuchsian differential equations of 
quite large orders, such as those found for the third and fourth contribution (x^''' and x'*'') to the 
magnetic susceptibility of square lattice Ising model. We use the previous connection matrices to get 
the exact explicit expressions of all the monodromy matrices of the Fuchsian differential equation 
for x'^' (and x'*') expressed in the same basis of solutions. These monodromy matrices are the 
generators of the differential Galois group of the Fuchsian differential equations for x^'^^ (and x^'^^)! 
whose analysis is just sketched here. 
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1 Introduction 

Since the work of T.T. Wu, B. M. McCoy, C.A. Tracy and E. Barouch [J, it is known that the expansion 
in n-particle contributions to the zero field susceptibility of the square lattice Ising model at temperature 
T can be written as an infinite sum: 

oo 

x(r) = ^x("n7^) (1-1) 

n=l 

of {n — l)-dimensional integrals |2 El El El El EI j the sum being restricted to odd (respectively even) n 
for the high (respectively low) temperature case. 

As far as regular singular points are concerned (physical or non-physical singularities in the complex 
plane), and besides the known s = ±1 and s — ±i singularities, B. Nickel showed 6 that x*-^"^^-* is 
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singular for the foUowing finite values of s = sh{2J/kT) lying on the \s = 1| unit circle (to = fc = 
excluded) : 

2-(s + -) = + i- + + — 

m2»+i =1, -n < m, k < 71 (1.2) 

In the following we will call these singularities: "Nickel singularities" . When n increases, the singularities 
of the higher-particle components of x(s) accumulate on the unit circle \s\ = 1. The existence of such 
a natural boundary for the total x(s), shows that x(s) is not D-finite (non holonomic^ as a function of 
.). 

A significant amount of work had been performed to generate isotropic series coefficients for x*^"^ 
(by B. Nickel 013 up to order 116, then to order 257 by A.J. Guttmann and W. Orrick^). More 
recently, W. Orrick et al. [S], have generated coefficients^ of x(s) up to order 323 and 646 for high 
and low temperature series in s, using some non-linear Painleve difference equations for the correlation 
functions |S1 1^1 El HH IT^ . As a consequence of this non-linear Painleve difference equation, and the 
remarkable associated quadratic double recursion on the correlation functions, the computer algorithm 
had a 0{N^) polynomial growth of the calculation of the series expansion instead of an exponential 
growth that one would expect at first sight. However, in such a non-linear^ non-holonomic, Painleve- 
oriented approach, one obtains results directly for the total susceptibility xi^) which do not satisfy any 
linear differential equation, and thus prevents the easily disentangling of the contributions of the various 
holonomic x'-"-''s. 

In contrast, we consider here, a strictly holonomic approach. This approach jl3l 1141 ITK] enabled us 
to get 490 coefficients^ of the series expansion of x^'^'' (resp. 390 coefficients for x^^"*); from which we 
have deduced |l3lll4[[T31ll6j the Fuchsian differential equation of order seven (resp.ien) satisfied by x^^-* 
(resp. x^^^)- We will focus, here, on the differential Galois group of these order seven and ten Fuchsian 
ODE'S. 



2 The Fuchsian differential equations satisfied by x^^\w) and 

Similarly to Nickel's papers [SI E], we start using the multiple integral form of the x'""'''Sj or more 
precisely of some normalized expressions x*'"'' ■ 

X^'^Hs) = ^± X^")(s), n = 3,4, ... (2.1) 

S+ = , T>Tc {n odd) 

s 

= (1-5-^)1/-^, T<Tc {n even) 

where: 

x(")(w;) = Jd"V (n^.j (2.2) 



\i=l 



with (each angle varying from to 2tt): 

n — 1 , , n 



d^v^ with y</). = o, i?H=i + nLii^ 



^The fact this natural boundary may be a "porous" natural frontier allowing some analytical continuation through it 
is not relevant here: one just need an infinite accumulation of singularities (not necessarily on a curve ...) to rule out the 
D-finite character of x- 

^A.J. Guttmann and W. Orrick private communication. 

*The short-distance terms were shown to have the form (T — ■ {log\T — Tc\)'' with p > ■ 

^We thank J. Dethridge for writing an optimized C-l — h program that confirmed the Fuchsian ODE we found for x^^^i 
providing hundred more coefficients all in agreement with our Fuchsian ODE. 
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^^"^^n^T-^-in^Pr^ (2.3) 

Instead of the usual[ni[7] variable s, we found it more suitable to use w — 55/(1 + s^) which has, by 
construction, Kramers- Wannier duality invariance (s <-> 1/s) and thus allows us to deal with both limits 
(high and low temperature, small and large s) on an equal footing ^13. .l^ #15) . The quantities Xj and 
yj can be written in the following form |13l 1141 ITS] : 

2w 



1 — 2w cos 0j + Y^(l — 2wcos(/)j)2 — ' 

2w , , 

yj = , (2.4) 
^ V(l-2wcos0j)2-4u;2 ^ ^ 

It is straightforward to see that x^"^ is only a function of the variable w. From now on, we thus focus on 
X^") seen as a function of the well-suited variable w instead of s [HIE]. One may expand the integrand 
in H2.2|l in this variable w, and integrate the angular part. 

We do not recall, here, the concepts, tricks and tools that have been necessary to generate very large 
series expansion for x^^H"^) X^^H"^) with a polynomial growth of the calculations jl3l 1141 IT^ . 

Given the expansion of x^^H"^) to w^^^, the next step amounts to encoding all the numbers in 
this long series into a linear differential equation. Note that such an equation should exist though its 
order is unknown^. Let us say that, using a dedicated program for searching'' for such a finite order 
linear differential equation with polynomial coefficients inw, we succeeded finally in finding the following 
linear differential equation of order seven satisfied by the 490 terms we have calculated for x'-^^: 

E «" • = ^ith: (2.5) 



dw 

a„ = w"-(l-4u;)^("-2)(^^4^)9(n-4)p^(^^^ 71 = 6, 5, • • • , 
where: 6{m) — sup{m, 0), and: ay = 

• (l-w) (l + 2w) (l-4u;)^(H-4w)^ {l + 3w + iw^) Priw) 

where Pt{w), Pq{w) ■ ■ •, Po{w) are polynomials of degree respectively 28, 34, 36, 38, 39, 40, 40 and 36 
in w |13j. 

Furthermore, besides the known singularities H1.2|) mentioned above, we remark the occurrence of the 
roots of the polynomial Py of degree 28 in w, and the two quadratic numbers roots of l + 3w + Aw'^ = 
which are not 17 Nickel singularities (they are not of the form (|1.2|l 'l. The two quadratic numbers 
are not on the s-unit circle : \s\ = \f2 and |s| = These quadratic numbers do not occur in 

the "physical solution" x^'^''- For ^7: near any of its roots, all the local solutions carry no logarithmic 
terms and are analytical since the exponents are all positive integers. The roots of P^ are thus apparent 
singularities |18lll9j of the Fuchsian equation H2.5|l . 

The order seven linear differential operator L7 associated with the differential equation satisfied by 
X^^-* has the following factorization properties |13l 1141 [T^ : 

L7 = © ie, ^6 = 1^3 • ^2 • (2.6) 

where* Li is a first order differential operator which has the first contribution to the magnetic suscep- 
tibility, namely x^^^ = 2w/(l — 4w), as solution. 

®A lower bound for the order of this linear differential equation would be extremely useful : such a lower bound does 
not exist at the present moment. 

'^Note that we, first, actually found an order twelve Fuchsian ODE and, then, we reduced it (by factorization of the 
differential operator) to a seventh order operator. This order twelve differential equation requires much less coefficients in 
the series expansion to be guessed than the order seven Fuchsian ODE we describe here ! ft is actually easier to find the 
order twelve differential equation than the order seven ODE !! 

SFor the notations see mO^Iig for x'^\ and for x^'*'- 
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In the same way, we found that the order ten linear differential operator Lio, associated with the 
differential equation satisfied by x^^\ has the following factorization properties |15l I16| : 

Lio = No(BLs, Ls = M2 ■ G(L) (2.7) 

where A^o is an order two differential operator which has the second contribution to the magnetic sus- 
ceptibility, x^"^^ as solution and where G{L) is an order four differential operator that can be factorized 
in a product of four order one differential operators ^Hl • 



3 Differential Galois group 

A fundamental concept to understand (the symmetries, the solutions of) these exact Fuchsian differential 
equations is the so-called differential Galois group UHl 1211 122L 1231 1^ , which requires the computation 
of all the monodromy matrices associated with each (non apparent) regular singular point, these ma- 
trices being considered in the same basis^. Differential Galois groups have been calculated for simple 
enough second order, or even third order, ODE's. However, finding the differential Galois group of such 
higher order Fuchsian differential equations (order seven for x^^h order ten for x^^^) with eight regular 
singular points (for x*"^') is not an easy task. Along this side a first step amounts to seeing that the 
corresponding (order seven, ten) differential operators do factorize in smaller order differential operators, 
as a consequence of some rational and algebraic solutions and other singled out solutions ^Hl- These 
factorizations yield a particular block-matrix form of the monodromy matrices jl6) . The calculation 
of local monodromy matrices in some "well-suited" local (Frobenius series solution) bases is easy to 
perform, however the calculation of the so-called connection matrices corresponding to the "matching" 
of the various well-suited local bases associated with the various regular singularities is a hard non-local 
problem. Of course from the knowledge of all these connection matrices one can immediately write the 
monodromy matrices in a unique basis of solutions|lf)j. 

From exact Fuchsian ODE's one can calculate very large series expansions for these (well-suited local 
Frobenius) solutions, sufficiently large that the evaluation of these series far away from any regular 
singularity can be performed^" with a very large accuracy (400, 800, 1000 digits ...). As far as x*^'^-' is 
concerned one can reduc e |TB] the calculation of these connection and monodromy matrices, to the 6x6 
matrices of an order six jl6| differential operator Lq appearing in the decomposition (|2.6|l . Connecting 
various sets of Frobenius series-solutions well-suited to the various sets of regular singular points amounts 
to solving a linear system of 36 unknowns (the entries of the connection matrix). We have obtained 
these entries in fioating point form with a very large number of digits (400, 800, 1000, ...). We have, 
then, been able to actually "recognize" these entries obtained in floating form with a large number of 
digits jl6j . 

In particular it is shown in |16j that the connection matrix between the singularity points and 1 /4 
(matching the well-suited local series-basis near w = and the well-suited local series-basis near w — 
1/4) is a matrix where the entries are expressions in terms of ^/3, tt, I/tt, I/tt^, ... and a (transcendental) 
constant introduced in equation (7.12) of 1 : 

1 r r r ... ... f vi-^ \"\y-^ = 



dy\ dy2dyi{^- 



2.- A A A ^''^"'^'''^{yl-i){yl-i)^ 

= .000814462565662504439391217128562721997861158118508 • • • 

Y ^ - 

(yi + 2/2) (2/2 + ys) {yi + y2 + ys) " 

This transcendental constant can actually be written in term of the Clausen function GI2 '■ 
- (y + 2 - 3V3 . Cl2q)) (3.1) 



^ These monodromy matrices are the generators of the monodromy group whi ch identifies with the differential Galois 
group when there are no irregular singularities, and, thus, no Stokes matrices )25| . 



^•^Within the radius of convergence of these series. 
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where CI2 denotes the Clausen fmiction 



cm = 

n=l 



sin(n ( 



This constant can also be written in terms of dilogarithms, polygamma functions or Barnes G- 
functions : 



4 



(1/6 



3\/3 



^ • Im {dilog (1/2- l/2i%/3)) 
• (1, 2/3) + VI/ (1, 5/6) - * (1, 1/6) - * (1, 1/3)) 



• (hi(2) +ln(7r) - 6 In 

2 TT 



G(7/6) \ 
G(5/6); 



The 6x6 connection matrix C(0, 1/4) for the order six differential operator Lq matching the Frobe- 
nius series-solutions around w = and the ones around w = 1/4, reads: 



C(0, 1/4) 



(3.2) 



29 
16 



1 














" 


1 





9\/3 
64Tr 














37r\/3 
32 














5 




3V3 
647r 








1 

167r2 


5 
4 


37r\/3 
32 


45 

2567r 





1 

32 





3 


157r%/3 
64 


225V3 37r\/3 
10247r 64 


64 









Not surprisingly-'^^ a lot of tt's "pop out" in the entries of these connection matrices. We will keep 
track of the tt's occurring in the entries of connection matrices through the introduction of the variable 
a — 2 in. 

The local monodromy matrices can easily be calculated |16) since they correspond, mostly, to "loga- 
rithmic monodromies" and will be deduced from simple calculations using the fact that each logarithm 
(or power of a logarithm) occurring in a (Frobenius series) solution, is simply changed as follows : ln(K;) 

ln(w) -|- r2, where fl will denote in the following 2 in. From the local monodromy matrix Loc{Q), 
expressed in the w = 1/4 well-suited local series-basis, and from the connection matrix H3.2|l . the 
monodromy matrix around w = 1/4, expressed in terms of the (w — 0)-well suited basis reads |16) : 



= C(0, 1/4) • Loc{^) ■ C(0, l/4)-i 



A 
B C 



(3.3) 



^^One can expect the entries of the connection matrices to be evaluations of (generalizations of) hypergeometric functions, 
or solutions of Fuchsian differential equations. 
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where 



A 
B 



and [ C ] read respectively: 



-24 a* 
-48 a" 


-48 pi 
(87 + 8a2)a4 





24 


32flp2 
4(75- 4^2) Q,2^ 







-300 a^n 
3(4a2 - 75)a^n 



with pi = 5 a* 



8 f72Q2^ p2 = 4 f7 - 75 - 15 and pg = 5 

24 Qf^ ~3Ma^n 1536 

24^4 -192 a^j^ 

24: a'^ 



AQ + ifla'^, and: 



C 



Note that the transcendental constant has disappeared in the final exact expression of H3.3|l which 
actually depends only on a and fl. This (a, J7) way of writing the monodromy matrix H3.3|) enables to 
get straightforwardly the A^— th power of (|3.3() : 



M^=o(l/4)(a, r!) 



N 



M„=o(l/4)(a, iV-f]) 



(3.4) 



Let us introduce the following choice of ordering of the eight singularities, namely oo, 1, 1/4, wi, 
-1/2, -1/4, 0, W2 {wi = (-3 +iy^{7))/8 and W2 — wl are the two quadratic number roots of 
1 + 3 w + 4ii;2 = 0), the first monodromy matrix Mi is, thus, the monodromy matrix M^=o{oo) (see 
(13.31) ) at infinity with a = f2 = 2i tt, M(oo), the second monodromy M2 matrix being the monodromy 
matrix at w — 1, A4(l), ... This is actually the particular choice of ordering of the eight singularities, 
such that a product of monodromy matrices is equal to the identity matrix^"^: 



Ml ■ M2 ■ M3 ■ Mi ■ Mi ■ Me ■ Mj ■ Ms = Id 
= M{oo)-M{l)-M{l/4)-M{wi) 
xM{-l/2) ■ 7W(-l/4) • M{0) ■M{w2) 



(3.5) 



It is important to note that relation H3.5|l is not verified by the (a, fl) extension (like H3.3|l ) of the 
monodromy matrices Mj. If one considers relation (|3.5() for the (a, fi) extensions of the Mi's, one will 
find that (|3.5|l is satisfied only when a is equal to il, but (of course^ ) this a = fl matrix identity is 
verified for any value of ft, not necessarily equal to 2i7r. 



3.1 Mutatis mutandis for x^^^ 

Similarly to x^'^-' the differential operator for x^^^ presents remarkable factorizations that yield a particu- 
lar block-matrix form of the monodromy matrices • Similarly, again, one can consider the (Frobenius 
series) solutions of the differential operator associated with x'^^^ around x — Aw"^ — and around 
the ferromagnetic (and antiferromagnetic) critical point x = 1 respectively. Again the corresponding 
connection matrix (matching the solutions around the singularity points a; = and the ones around the 
singularity point x — — ^) have entries which are expressions in terms of tt^, rational numbers 

but also of constants like constant 
Riemann zeta function, as follows : 



introduced in ^ which can actually be written in term of the 



_l£2r3 



'47r2 



C(3)) 



(3.6) 



^^Of course, from this relation, one also has seven other relations deduced by cyclical permutations. 

^^A matrix identity like 13.51 yields a set of polynomial (with integer coefficients) relations on H = 2i7r. The number 
TT being transcendental it is not the solution of a polynomial with integer coefficients. These polynomial relations have, 
thus, to be polynomial identities valid for any Q. 



6 



The derivation of the two results Ij^-lfl . H3.6|l for the two transcendental constants and has never 
been published^'' but these results appeared in a conference proceedings [201 • We have actually checked 
that and we got in our calculations of connection matrices displayed in \27\ as floating numbers 
with respectively 421 digits and 431 digits accuracy, are actually in agreement with the previous two 
formula. These two results (|3.1(l , (|3.6() provide a clear answer to the question of how "complicated and 
transcendental" some of the constants occurring in the entries of the connection matrices can be. These 
two remarkable exact formulas (|3.1|1 . H3.6|l are not totally surprising when one recalls the deep link between 
zeta functions, polylogarithms and hypergeometric series 28, 29, 30, 31 . Along this line, and keeping in 
mind that we see all our Ising susceptibility calculations as a "laboratory" for other more general problems 
(Feynman diagrams, ...), we should also recall the various papers of D. J. Broadhurst | 32ll33[l5Hl35| where 
C^2(f ) and C(3) actually occur in a Feynman-diagram-hypergeometric-polylogarithm-zeta framework 
(see for instance equation (163) in |32j ) . 

Similarly to the previous results for x^'^^ the monodromy matrices written in the same basis of 
solution, deduced from the connection matrices and the local monodromy matrices are such that a 
product in a certain order of them is the identity matrix. Denoting by Mx=o{0): Mj;=o(l), Me=o(4) and 
Mx=o{oo) the monodromy matrices expressed in the same a: = well-suited basis, one obtains: 

M,=o(oo) • A/,=o(4) • M^=o{l) ■ Mx=o(0) = Id (3.7) 

This matrix identity is valid irrespective of the "not yet guessed" constants |16| . 



4 Conclusion 

The high order Fuchsian equations we have sketched here present many interesting mathematical prop- 
erties close to the ones of the so-called rigid local systems j36| . these rigid local systems exhibiting 
remarkable geometrical interpretations |37) as periods of some algebraic varieties. This "rigidity^'^" 
emerges through the log-singularities of the solutions of these Fuchsian ODE's: the powers of the log- 
arithms of these solutions are "smaller" than one could expect at first sight. It is worth noting that 
almost all these mathematical structures, or singled-out properties, we sketched here, or in previous 
publications p2lll4' 15, 16 , are far from being specific of the two-dimensional Ising model : they also 
occur on many problems of lattice statistical mechanics or, even^®, as A. J. Guttmann and I. Jensen 
saw it recently, on enumerative combinatorics problems like, for instance, the generating function of the 
three-choice polygon (38) . 

We have also seen in some of our calculations ^^E] a clear occurence of hypergeometric functions, 
hypergeometric series and in some of our calculations (not displayed here) generalizations of hyperge- 
ometric functions to several complex variables: Appel functions 39 , Kampe de Feriet, Lauricella-like 
functions, polylogarithms 32, 33, 34, 35 , Riemann zeta functions, multiple zeta values, ... The occurence 
of Riemann zeta function or dilogarithms in the two remarkable exact formulas (|3.1|l . H3.6|l is not to- 
tally surprising when one recalls the deep link between zeta functions, polylogarithms and hypergeometric 
series gHl E3 EOl EH ■ 

We think that such "collisions" of concepts and structures of different domains of mathematics 
(differential geometry, number theory, ...) are not a consequence of the free-fermion character of the 
Ising model, and that similar "convergence" should also be encountered on more complicated Yang- 
Baxter integrable models^'', the two-dimensional Ising model first "popping out" as a consequence of 
its simplicity. In a specific differential framework some of these interesting mathematical properties can 
clearly be seen in the analyzis of the differential Galois group of these Fuchsian equations. 

We have underlined the fact that, beyond a general analyzis of the differential Galois group one 
can actually find the exact expressions of the non-local connection matrices from very simple matching of 
series calculations, and deduce, even for such high order Fuchsian ODE's, explicit representations of all 

^*We thank C. A. Tracy for pointing out the existence of these two results 13.11 . 13. 6i and reference 1261 . 
^^Let us recall that hypergeometric functions are totally rigid. 

^®The wronskian of the corresponding differential equation in 1381 is also rational, the associated differential operator 
factorizes in a way totally similar to the Fuchsian ODE's for x^^^ a-^d x^^' i large polynomial corresponding to apparent 
singularities also occur, ... 

'^^The comparison of the Riemann zeta functions equations obtained for the XXX quantum spin chain 1401 with the 
evaluations of central binomial in [32j provides a strong indication in favor of similar structures on non-free-fermion Yang- 
Baxter integrable models. 
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the monodromy matrices in the same (non-local) basis of solutions, providing an effective way of writing 
explicit representations of all the elements of the monodromy group. The remarkable form, structures 
and properties (see (|3.2|) . H3.3|) . H3.4|l . (|3.5|l 'l of the monodromy matrices in the same (non-local) basis of 
solutions is something one could not suspect at first sight from the general description of the differential 
Galois group. 
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